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I) : 

^ , Abstract 

Open string field theory in the level truncation approximation is considered. It 
V^ . is shown that the energy conservation law determines existence of rolling tachyon 

lO ! solution. The coupling of the open string field theory action to a Friedmann- 

p^ ' Robertson- Walker metric is considered which leads to a new time dependent rolling 

^ ■ tachyon solution is presented and possible cosmological consequences are discussed 

o 

00 : 1 Introduction 

o 

Consideration of fundamental theories such as M/String Theory in the cosmological con- 
^ ■ text continues to attract attention in the literature. One of the interesting questions is the 

role of tachyon in String Theory and Cosmology. The great progress in our understand- 
ing of tachyon condensation was made in the past decade [1-3], but a lot of interesting 
issues are still open. Among the most important ones is a better understanding of the 
dynamics of tachyon condensation process. 

In this context many works have been devoted to the study of time-dependent so- 
lutions. Probably one of the most fascinating frameworks for this is Open String Field 
Theory (OSFT) [4] which reattracted a lot of attention after a recent work [5]. Despite 
the recent renewal of interest to OSFT no smooth solutions interpolating between two 
inequivalent vacua even at the lowest level truncation order [6, 7] were found. One of 
the exceptional features of the level truncation approach is that corresponding action 
contains infinitely many time derivatives, i.e. it is non-local. Resulting models have a 
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rich set of properties that might be essential for the development of stringy cosmology. 
Substantial investigation of this type of models was performed [8,9, 18-33,35-37,42,43] 
after seminal paper [8] where many useful issues were revealed, among which the presence 
of rolling solution with widely increasing oscillations. Further research of time-dependent 
rolling solutions was also performed in particular in [10]. Recently several time dependent 
rolling tachyon solutions in OSFT were found [12,13], which confirmed puzzling behavior 
of the solutions found earlier [8,10,11]. In these works investigations were performed 
in usual space-time coordinates. At the same time is it known that in the light-cone 
gauge the theory becomes local in light-cone time [14]. In fact the question of the iden- 
tity of light-like and time-like cases is under investigation. In the light-like case (with 
dilatonic damping), the gradient flow forces tachyon to asymptote to the true vacuum at 
late times [15]. In the time-like case, rolling configuration between two non-equivalent 
vacua is forbidden by energy conservation law in Minkowski space time if one considers 
only tachyon excitation in the level truncation approximation [8,16]. Investigations in 
the last case though did not take into account any effects of gravity, this makes such an 
investigation inconsistent from cosmological point of view. 

It turns out that if we consider the same action coupled to the gravity in Friedmann- 
Robert son- Walker (FRW) metric the situation changes: there appears a tachyon solution 
which tends to the true vacuum at late times [16]. It is interesting to note that because 
dilaton appears from the same string sector as graviton, inclusion of dilaton into the 
tachyon action can qualitatively reproduce behavior of the tachyon in curved space. Thus 
obtained solution has accordance with the results obtained in light-like and time-like cases 
at least at the lowest level approximation. 

In the present work we consider scalar field dynamics with infinitely many time deriva- 
tives minimally coupled to the Minskowski and Friedmann-Robertson- Walker gravita- 
tional backgrounds. 

The structure of the work is the following. In the first section we will give a brief 
introduction and physical motivation. In the second section the model which appears 
from Open String Field Theory will be presented and problem of existence of interpolating 
solutions between two inequivalent vacua will be discussed. In the third section we will 
consider the model minimally coupled to gravity and demonstrate an intriguing difference 
compared to the case without gravity: the existence of desired solution. Numerical 
techniques will be described in Sec. 4. Results of numerical calculations will be presented 
in Sec. 5. Finally we will summarize main results. 

2 The Model 

The action of bosonic cubic string field theory has the form 



1() J '^ >j 
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where Qq is the open string coupling constant, Qb is BRST operator, * is noncommutative 
product and $ is the open string field containing component fields which correspond to 



all the states in string Fock space. 

Considering only tachyon field (l){x) at the level (0,0) the action ([2]) becomes 
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where ot is the Regge slope, K = ^j^, is a scalar field, $ = e "0, k = a'lnK, 
n = -^^d^^J—gg^^dy and A = \K~^ was added to the potential to set the local minimum 
of the potential to zero according Sen's conjecture [38]. In what follows we will work in 
units where g^ = 1. 

The action ([2]) leads to equation of motion 



The Stress Tensor for our system i^j 
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+2k I dp {d^e^P^K^<^'^){dpe~^P^^). 
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The energy is defined as E{t) = T^^ and pressure as P{t)i = — T/ (no summation) and 
for our system are 



S — Sk + Sp + A + Snll + Snl2, V — Sk — Sp — A — Snll + S., 



nl2 



where 



s, = j{dc^r, Sp = -\<f>' + ^^^ 

S^i^ = k [ dp {e''P°K'^<^^) {-De'''P°^) , 
Jo 

Sni2 = -k [ dp {de^P^K"^^^) {de~^f'°^) . 

In this paper we will be interested in spatially homogeneous configurations for which 
Beltrami-Laplace operator used above takes the form D^ = —d^ and nonlocal operator 
e'^'^" becomes e'^°^ = e~^^ . 



^Note that here and below integration over p is understood as hniit of the following regularization 



dpfip) = hni Hni 

Q £1^ + 62^ + 



dpfip). 



The symbol e^^ Lp comprehend aqj 

e''''^{t) = C,mt) (5) 

where 

Nonlocal terms Eni^ and £ni2 contain e~''^^ which as it can be easily seen might lead to 
the growing kernel in the integral representation of the non-local operator, that is why 
we will try to avoid calculation of e"*^^^ and will use the following representation for 
nonlocal energy terms Snii and Sni2 which are valid on the equation of motion for the 



scalar field 



Sni2 = -k f dp (d{-a'd^ + l)e(2-'')'=^'$) fde'^P^'^ 



2.1 Energy conservation 

Energy conservation laws always have the deep physical sense, note that because of pres- 
ence of infinitely many time derivatives we need to prove energy conservation explicitly. 
Similar investigation was performed in [8], although in that work it was used represen- 
tation of pseudo-differential operator via the summation over infinite series expansion 
with which one always needs to be very careful with regard to convergence issues. The 
approach used here is advantageous from the point of view of numerical calculations, 
because in order to define action of the exponential operator we need to do only one well 
defined integration. 

Taking into account that models with different types of potentials are currently under 
the consideration in the literature [29,31,43] we will show the energy conservation for 
arbitrary potential. 

Claim 1. u The Energy 

/ -I rl 

E = —{d(j)f - -02 + y($) +A + k dp {{-a'd^ + l)e(2-p)'=9'$)'9'(9e'=''^'$), 
2 2 Jq 

is conserved on the solutions of equation of motion 

.2 , ,,.w^^_dV{<l>) 
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^It is easier to use integral representation for the operator e^ while considering Minkowski back- 
ground, in the FRW case though it becomes impossible to generalize such an approach and we define 
the operator e^ in terms of solution of the boundary value problem for diffusion equation, see [31] for 
the details. 

•^For simplicity we will use symbolic notation for nonlocal operator e keeping in mind that it is in fact 
defined by ([5]) , also integration over p must be understood as limit of the corresponding regularization 
as indicated earlier. 



of the corresponding actions 
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where V{^) is any polynomial potential and A d B = AdB — BOA. 
Proof. We compute 
dE{t) 



la^i^'i. 



dt 
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Using following identity (for details on its derivation see Appendix) 



dp{eP^"ip)d^e^'-P^^"(j)) = ipe^"(j), 



equation of motion and definition of field $, we have 

^^ = a'd(j)d''<P - 4>d(t> + ^5$ + d^7^\a'd^ - l)e'=^'$ 
dt ^ ^ ^ ^ 5$ 
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In the next paragraph we consider one important physical consequence of this energy 
conservation law. 



2.2 Existence of the Rolling Solution 

We already indicated that from physical perspective we are interested in solutions inter- 
polating between two inequivalent vacua. We start our consideration by looking for sta- 
tionary configurations $o- Substituting it into equation of motion ^ we get $o = -^^"^O' 
which has two constant solutions: $o = and $o = K~^ . We should thus be looking 
for solutions interpolating between those stationary points. The following claim though 
shows that energy conservation forbids existence of such solutions. 

Claim 2. □ There do not exist continuous solutions of equation ^ which satisfy boundary 

conditions 

^ / N I 0, t -^ oo, 

lim $(t) = { l_^ ^ (6) 



K- 



-oo 



or vice-versa (in terms t — ^ —t). 

"* Similar claim for p-adic string model was proved in [8], which rules out the possibility that tachyon 
may roll monotonically down from one extremum reaching the tachyon vacuum. 



Proof. Let us assume existence of such solution and calculate energy at the extremum 
points, we get £'($ = 0) = A and £'($ = K~^) = —^K~^ + A, i.e. energy values at 
t — i> +00 and t -^ —00 are different what due to conservation law rules out existence of 
solutions satisfying (jS]). D 

As we can see energy conservation law plays crucial role in the existence of the time 
dependent solutions in the level truncation approximation to OSFT. The above statement 
could potentially be generalized to the case of full OSFT because for the action with 
cubic interaction solution interpolating between maximum and minimum in the effective 
potential has to interpolate between vacua with different energy. 

3 The Model Coupled to the Gravity 

In this section we would like to consider tachyon dynamics in FRW background what 
allows us to take into account gravity effects and makes the research more consistent 
from cosmological point of view. Consider the model 

S = ^J d'xV^g {^R + \ct.U,4> + \4>' - \K'e - a) , (7) 

here m? = g^M^i and we will work in units where a' = 1. As a particular metric we will 
consider the FRW 

ds'^ = —dt^ + a^{t){dx\ + dx\ + dxl), 

for which the Beltrami-Laplace operator for spatially-homogeneous configurations takes 
the form D^ = —d"^ — 3H{t)d = —2)^. Scalar field and Friedmann equations are 

3H^ = ^S, 3H^ + 2H= -^V. (9) 

Inclusion of the gravity considerably modifies the dynamics of the system. In Minkowski 
background we have shown that energy conservation law forbids dynamical interpolation 
between two inequivalent vacua. In FRW background energy of the scalar field alone 
does not conserve any more due to the Hubble term and as a result the restrictions on 
existence of such solutions are removed. 

To find boundary conditions for possible solutions let us consider constant scalar field 
solution, $0- In this case the scalar field equation (jHj) becomes 

$0 = K'^<!>1, (10) 

and first equation in (Q 

3H', = ^S{<l>o). (11) 
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Figure 1: The potential 

Equation (fTOl) has two solutions: $0i = and $02 = K~^, substituting which into fITT]) 
we obtain corresponding values for a Hubble parameter Hq-^ = (ISK^)'^^"^ and Hq^ = 0. 
Note that from cosmological perspective we are interested only in positive values for 
the Hubble function, so we can expect rolling solutions with the following boundary 
conditions 



lim$(t) 




lim H(t) 




"1/2 



t -^ 00, 
t — >• —00, 



;i2) 



or vice- versa (in terms of t — > —t). 

To analyze physical situation let us consider potential in which motion is expected. 
Naive extraction of potential from the model action ([7]) results in V{^) = — 1$^ + 
|i^'^$^ -|- A. The constant A represents the D-brane tension which according to Sen's 
conjecture must be added to cancel negative energy appearing due to the presence of 
tachyon. We obtained two types of solutions. The first one is an ordinary rolling solution 
which starts from $ = and goes towards configuration $ = K~^ which is associated 
with the true vacuum. This solution can be interpreted as a description of the D-brane 
decay. The second one is a rolling solution which goes in the opposite direction, which 
appears in this model possibly because of the non-locality in the interaction. It is known 
that nonlocal dynamics has many interesting properties which are not possible in the 
local case. In particular the "slop effect" [8,21,27] which we can observe in the obtained 
solutions (Fig|2l [3]) when the scalar field goes beyond the values from which the scalar 
field configuration starts - situation which is not possible in the local models. Potentially 
a similar effect can initiate non-symmetry in the potential in ekpyrotic [39] and cyclic 
cosmology [40]. 

4 Numerical Solution Construction 

In order to construct numerical solution we operate with scalar field equation of motion 
([8]) and the difference of equations Qj, specifically we solve 



-Vl + l)e2^^H$ 



3^2 



K-'<l> 



'13a) 



H 



2m2 



{V + £). 
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We discretize equations fll3ap - fll3bp by introducing a lattice in t variable and then solving 
resulting system of nonlinear equations using iterative relaxation solver controlling error 
tolerance with discrete L2 and Loo norms. 



4.1 Discretization 



When solving discretized equations fll3al) - fll3bl) the nontrivial part from computational 
point of view is efficient evaluation of e^'^^^Hcl) for p E [0,2]. This operator could be 
interpreted in terms of initial value problem for the following diffusion equation with 
boundary conditions [31] 



dpip{t, p) = dtip{t, p) + 3H{t)dtip{t, p), 



(14) 



Once solution of this equation is constructed we have e^'^''^H$(t) = Lpi^p^t). 

To solve flT^ we used second order Crank-Nicholson scheme which is based on the 
following stepping procedure 



V9(t,p + A^ 



1 + kApVl 



1 - k\Vl 



-1 



^ + o(Aj| 



1)^1 



where D^ denotes discretization of D^ (it thus has a finite norm) and Ap is a step size 
along p variable. We used the following tri-diagonal discretization scheme for Dy 
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where A^ is the step size along t and discretization of Hubble function is H^ = H{kAt). 
This is a usual symmetric discretization scheme on the uniform lattice modified on the 
interface to guarantee boundary condition fll3bp for smooth solutions which tend to 
constants as t — *> ±00. 



4.2 Comparison the different methods 

In order to exclude possible artifacts of the specific numerical scheme described above we 
tried Chebyshev-pseudospectral method which is known to generally have exponential 
convergence [41]. Such scheme is known to have very different properties [41] compared 
to finite difference scheme described above, but it produced the same results up to the 
approximation error which provides confidence in the existence of the rolling solutions 
reported in this work. 
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Figure 2: Solutions of the scalar field 
right) for mi = 1. 
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([8]) and Friedmann equation ([9]) $ and H (left to 



5 Rolling Tachyon Solution 

Solutions of ([8]) and ([9]) are presented on FigJJl [3l It is interesting to note that on 
these figures shapes of scalar field look very similar up to reflection over vertical axis 
while shapes of Hubble function are different. In order to ensure that this is not an 
artefact or error in numerical calculations we performed the following test. First, we 
performed reflection, t -^ —t, of the scalar field from Figl2] and compared in with scalar 
field from Fig|2l while difference was small it was well above the error tolerance. We also 
numerically computed discrepancies between left and right sides of equations of motion 
where Hubble function was from the solution while scalar field was taken from another 
solution and mirrored. While the discrepancy for scalar field was small, the discrepancy 
for Hubble function was very large, larger than the Hubble function itself. This test 
makes us confident that qualitatively different shapes of Hubble function presented on 
Fig 121 [31 are correct. 

Another interesting property of the solution on FiglSlis that while scalar filed generally 
rolls down at first it climbs up, we already noted this "slop effect" in the end of section 
131 In fact similar behavior is typical for solutions of nonlocal equations. In Minkowski 
metric such behavior was noted by many authors who used different numerical techniques 
[21,27,42]. It is interesting that in FRW metric this property is preserved. 
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Figure 3: Solutions of the scalar field 
right) for mi = 1. 
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Figure 4: Solutions of the scalar field ([H]) and Friedmann equation ([9]) $, H ioi rrip 
0.1,1,10 (left to right). 



5.1 Solutions for different rup 

Motivated by the fact that string scale does not exactly coincide with Planck mass and 
as a consequence there is some freedom in settling mi we investigated how the shape 
of solutions behave for different rrip. Solutions of ([8]), ([9]) for different values of nip are 
presented on FigJH and O 

We can see that the profiles for different m^ for the first case are very similar. It is 
interesting to note how on Fig l5] with the growth of rrii oscillation in the profile of Hubble 
function disappears. 



5.2 Particular case of j9-adic action 

The p-adic string model represents a popular toy-model that was proposed in [34]. The 
model contains one scalar field and nonlocal interaction. Formulated in arbitrary space- 
time dimensionality the model has a parameter p, which is initially taken to be a prime 
number. It was shown later the model describes quite well some physical processes, 
despite its intrinsic limitations. Recently this model started to attract attention as a 
cosmological toy model for describing inflation [35,36,43]. The jo-adic string model was 
later considered as being minimally coupled to gravity in Friedmann-Robertson- Walker 
metric. In this subsection we would like to show some intriguing properties within this 
context. 

The p-adic action is given by 
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The infinite number of space-time derivatives are manifest in the pseudo-differential op- 



erator p 2°, where □ = —d"^ + V^, p 
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Figure 5: Solutions of the scalar field ([H]) and Friedmann equation ([9]) $, H ioi nip 
0.1,1,10 (left to right). 

Considering action (fT5l) for p = 2 after the field redefinition tp = e~i^^^°(f) we have 



S, 



d'^xL2 



4 / d^'x 
9p 



4^^ + ^(e^°V^f 



(16) 



which looks rather similar to the tachyon OSFT action in the level truncation approxi- 
mation without kinetic term with only other difference in common signs. 
Indeed, if we denote the approximate lagrangian for the OSFT case 



^OSFTap 



V_1^3. /c□,^3 



we see that for the case of p-adic string model with p = 2 if we neglect the difference in 
the factors in exponential operator we have 



^p=2 



-L 



OSFTa^ 



This observation does not affect dynamic in the case of Minkowski space-time because 
sign in front of the Lagrangian does not enter scalar field equation. The situation changes 
crucially in Friedmann-Robertson- Walker background because the sign affects equation 
for Hubble. Thus if in the case of OSFT action in approximation neglecting kinetic term 
we had monotonically increasing Hubble function in the case of p-adic model for p = 2 we 
obtain monotonically decreasing Hubble function for the same scalar field configuration. 
Corresponding solutions are presented on FigjGl 

On of the issues is that the p-adic sting model in the FRW case is not currently 
known and the choosing the right sign is nontrivial, especially due to the absence of 
usual canonical kinetic term. 



11 





Figure 6: Solutions of the Friedmann equations $, H: negative sign in front of the 
Lagrangian (left), positive sign in front of the Lagrangian (right). 
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Figure 7: Energy (left) and pressure (right) for solutions from Fig|2l 

5.3 Energy and Pressure for Rolling Tachyon Solutions 

Usually when considering dynamics in nontrivial background such as FRW it is instructive 
to consider simpler case of Minkowski metric. The intriguing fact though is that there is 
no such a possibility since energy conservation law forbids in Minkowksi case solutions 
like were presented in Figj2l [31 Nevertheless it is possible to qualitatively characterize 
new features of the system. 

Dynamics of energy and pressure for the scalar field from Figl2] is presented on FigJTl 
We can see that dynamics is different from what one might expect in Minkowski case. 
We obtained solution for which both energy and pressure nontrivially tend to zero at 
large times. For the solution from Fig|3]the energy and pressure dynamics is different, 
see FigJHl The pressure starts from zero and goes to a negative constant, while the energy 
starts from zero and goes to a constant of the same value with positive sign. 

In both cases we obtained nontrivial partially negative pressure and like in [3] non- 
trivial equation of state what might be especially interesting in the light of cosmological 
applications. 

5.4 Cosmological Applications 

To discuss cosmological consequences it is interesting to refer to the concluding remarks 
and open questions of the work [8] which in fact initiated a whole series of investigation of 
nonlocal dynamics in the models with infinitely many time derivatives [8]. Original paper 
considered string field theory and p-adic theory in Minkowski space time. It was noted 
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Figure 8: Energy (left) and pressure (right) for solutions from FigJ31 

that the most puzzling result from physical point of view is that none of the solutions 
obtained there appeared to represent tachyon matter. In other words there were not 
found solutions where with varying values of energy the pressure tends to zero for long 
times. It is notable that solutions obtained in this work have such behavior. As we can 
see on FigJT] the pressure goes to zero for long times at the end of the evolution. In 
fact this represents first alternative of possible cosmological applications of the obtained 
solutions. 

The second alternative is to consider obtained solutions in the context of description 
of very early Universe. This investigation is motivated by the shape of the solutions with 
decreasing positive Hubble function presented at the FigJ21 The research in this direction 
was started in [35] based on p-adic string model and followed by [43] in which CSSFT 
string model was considered as well. 

The third alternative is to consider configuration from Figj3] in the context of late 
time acceleration. 
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6 Summary 

In this paper we have studied tachyon dynamics in string theory in Friedmann- Robertson- 
Walker background. The main results are: 

• It was shown that filed configurations which interpolate between vacua with dif- 
ferent scalar field energies are possible in the FRW background and corresponding 
rolling solutions were obtained. 

• We found nontrivial dynamics in the FRW background which differs from dynamics 
in Minkowski case. Particularly, we found interpolating solution which goes from 
maximum to the true minimum for which the pressure approaches zero at long 
times while energy varies starting from the constant which is equal to the brane 
tension and tends to zero at the end of the evolution. 

Among main results we have also considered in details the properties of constructed 
solutions and have shown that the evolution in tachyon potential is possible in both 
directions, and what is even more puzzling profiles for Hubble functions are different in 
those cases. In one direction Hubble function is an increasing almost monotonic function, 
while in another we observe significant oscillations during evolution. Moreover it was 
shown that by varying string scale we can change the shape of Hubble function and for 
some string/Planck scale ratio the oscillations disappear. We also considered dynamics 
in p-adic string model for particular value of a parameter p = 2 and have shown that 
dynamics in the FRW background drastically differs from the corresponding one in usual 
string theory. 

It is interesting to note that for numerical construction of solutions presented here 
we had to abandon explicit iterative techniques which were very successful in previous 
investigations of rolling tachyon solutions, see [31,43] and references therein. Instead 
we used slower but more generic relaxation methods with Crank-Nicholson scheme to 
compute nonlocal operators. 
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Appendix 

In this appendix we will prove the following identity 



dp(e^^V) d^ (e(i-^)^V) = (/^e^V, (17) 



where symbol e''^ ip comprehend as [16] 

It is well known, that for functions $(if:) which are continuous and bounded on the 
real axis the following identity have a place 

\im^C,mt) = ^t) (18) 

We can formulate the following lemma 



Lemma 1. For continuous and hounded functions ■il){t) and Lp{t) the following identity 

has a place 

1 

d^ 





where left side we understand as 



1 

j dpCp[^]{t)?C,.,mt) = ^C[ip, d^ = 



(19) 



lim lim / dpdpCp[^]{t)d'Ci^p[ilj]{t) 
and right hand side we understand as 

^ciij = p{t)c,mt) - c,[p]{t)m 

Proof. We have 

1 — £2 1 — £2 1 — £2 

32^ r;.n / J„/n2/ 



dpCMd'c^.pitP] / dpCMid'c^-pm- / dp{d'cM)'^C- 



a^p[ 



£1 



We will use the fact that for p > the function Cp[(p\{t) is a solution of the diffusion 
equation, i.e 

^Ci_,M(t) = -^Ci_,M(t), p>o 

The proofs of these identities follow from integral representation. 
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Taking into account the identities written above we get 

l-e2 1-62 l-e2 

ei £i £i 

1-52 



/ dp (c,[^]§-^c^^,m + ^C,[v.]Ci_,M 



£1 
1-62 



dp-icMc,^pm = -c^s^ivn^m + CsA^]Ci-eAi^] 

and flTSj) . we can take the hmit £i -^ +0, £2 -^ +0 and get flT9|) . D 

References 

[I] A. Sen, Int.J.Mod.Phys. A 20, 5513 (2005). 

[2] A. Sen, J. High Energy Phys. 0204, 048 (2002). 

[3] G. W. Gibbons, Phys. Lett. B 537, 1 (2002); G. W. Gibbons, Class. Quant. Grav. 
20, S321 (2003). 

[4] E. Witten, Nucl. Phys. B, 268, 253 (1986). 

[5] M. Schnabl, Adv. Theor. Math. Phys. 10, 433 (2006) [arXiv:hep-th/0511 286]. 

[6] V.A. Kostelecky and S.Samuel, Phys. Lett. B 207, 169; Nucl. Phys. B 336, 263 
(1990). 

[7] N. Moeller, A. Sen, B. Zwiebach, J. High Energy Phys. 08 (2000) 039. 

[8] N. Moeller and B. Zwiebach, J. High Energy Phys. 10, 034 (2002). 

[9] H. Yang, JHEP 0211 (2002) 007. 

[10] M. Fujita and H. Hata, Time dependent solution in cubic string field theory, JHEP 
0305, 043 (2003) |!arXiv:hep-th/0304163] . 

[II] M. Fujita and H. Hata, Phys. Rev. D 70, 086010 (2004) |arXiv:hep-th/0403031| . 



[12] M. Schnabl, Phys. Lett. B 654, 194 (2007) [arXiv:hep-th/0701248; 

[1 3] M. Kiermaier, Y. O kawa, L. Rastelh and B. Zwiebach, JHEP 0801, 028 (2008) 
|arXiv:hep-th/0701249| . 



16 



[14] M. Kaku and K. Kikkawa, Phys. Rev. D 10, 1110, 1823 (1974). 

[15] S. Hellerman, M. Schnabl, arxiv: 0803.1 184 [hep-th]. 

[16] L. Joukovskaya, arxiv 0803.3484[hep-th]. 

[17] V. Forini, G. Grignani, G. Nardelli, J. High Energy Phys. 0604:053, 2006. 

[18] I.Ya. Aref 'eva, L.V. Joukovskaya and A.S. Koshelev, J. High Energy Phys. 09 (2003) 
012. 

[19] Ya. Volovich, J. Phys. A 36, 8685 (2003). 

[20] L. Joukovskaya, Ya. Volovich, jmath-ph/0308034 1 

[21] V. Forini, G. Grignani, G. Nardelh, J. High Energy Phys. 0503, 079 (2005). 

[22] L.V. Joukovskaya, Theor. Math. Phys. (EngL TransL) 146, 335 (2006), 
larXiv:0708.0642[ 

[23] Vladimirov and Ya.I. Volovich, Theor. Math. Phys. (Engl. Transl.) 138, 297 (2004). 

[24] V.S. Vladimirov, Theor. Math. Phys. (Engl. Transl.) 149, 1604 (2006). 



[25] D.V. Prokhorenko, |arXiv:math-ph/0611068 



N. Barnaby, N. Kamran, JHEP 0802:008, 2008. 
[27] I.Ya. Aref eva, L.V. Joukovskaya, JHEP 0510 (2005) 087. 
[28] I.Ya Aref eva, AIP Conf Proc.826:301-311,2006. 
[29] G. Calcagni, J. High Energy Phys. 05 (2006) 012. 

[30] I.Ya Aref eva, L.V. Joukovskaya, S.V. Vernov, J. High Energy Phys. 07 (2007) 087. 
[31] L. Joukovskaya, Phys. Rev. D 76, 105007 (2007). 
[32] L. Joukovskaya, AIP ConfProc. 957, 325 (2007). 

[33] G. Calcagni, M. Montobbio, G. Nardelh, Phys.Rev.D 76, 126001 (2007). 
[34] L. Brekke, P.G. Freund, M. Olson and E. Witten, Nucl. Phys. B 302, 365 (1988). 
[35] N. Barnaby, T. Biswas, J. M. Cline, J. High Energy Phys. 04 (2007) 056. 
[36] J. E. Lidsey, Phys. Rev. D 76, 043511 (2007). 

[37] I.Ya. Aref eva, L.V. Joukovskaya, S.Yu. Vernov, arXiv:0711.1364[ hep-th]. 
[38] A. Sen, Int. J. Mod. Phys A14, 4061 (1999). 

17 



[39] J. Khoury, B.A. Ovrut, P.J. Steinhardt and N. Turok, Phys. Rev. D 64, 123522 
(2001). 

[40] P.J. Steinhardt, N. Turok, Science 296, 1436 (2002). 

[41] B. Fornberg, Cambridge University Press, 1996, ISBN 0521645646. 

[42] N. Moeller, arxiv: 0804.0697 [hep-th]. 

[43] D. Mulryne, N. Nunes, arxiv: 0805.0449. 



18 



